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ABSTRACT 


The  results  of  measurements  of  thermoacoustically  generated 
temperature  gradients  in  short,  thin  plates  located  in  a  resonant  tube 
are  reported.  The  temperature  gradient  results  from  a  heat  flow 
generated  by  the  acoustic  field.  Measurements,  at  both  axial  and 
transverse  positions  in  the  tube,  are  made  for  pressure  amplitudes 
ranging  from  approximately  145  dB  to  162  dB  re  20  pPa  (rms).  in 
argon  and  helium,  for  single  plates  and  stacks  of  up  to  five  plates 
having  separations  from  approximately  4  to  40  thermal  penetration 
depths.  The  results  are  compared  to  predictions  based  upon  a  theory 
by  Wheatley  and  others  (J.  Wheatley,  et  al.,  Journal  of  the  Acoustical 
Society  of  America,  Vol.  74,  pp.  153-170,  1983).  For  pressure 
amplitudes  below  150  dB,  the  measurements  agree  well  with  theory. 
At  higher  pressure  amplitudes,  the  agreement  diminishes.  It  is 
concluded  that  non-thermal  effects  are  responsible  for  the 
i  discrepancies. 
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Thermoacoustic  heat  transport  has  been  studied  in  detail  by  the 
late  John  Wheatley  and  associates  [Refs.  1  and  2].  One  instrument 
developed  in  these  studies  is  the  thermoacoustic  couple,  or  TAC. 
Wheatley  found  that,  when  a  TAC  is  placed  in  a  resonant  tube,  driven 
so  that  the  acoustic  wavelength  is  much  greater  than  the  length  of  the 
TAC,  a  temperature  gradient  is  induced  across  it.  The  gradient  is  a 
function  of  the  acoustic  pressure  amplitude,  the  type  of  gas  in  the 
tube,  the  thermodynamic  properties  of  the  TAC  plate  material,  the 
length  of  the  plate,  and  the  plate  configuration  (i.e.,  single  or  multiple 
plates). 

It  is  interesting  that,  while  Wheatley  states  that  the  thermoacous¬ 
tic  phenomena  can  be  observed  in  a  single  plate  [Ref.  l:p.  155],  his 
single  plate  results  indicate  a  large  deviation  from  theory  [Ref.  l:p. 
158].  This  discrepancy  is  noted,  but  there  is  no  report  of  further 
investigation  into  its  cause. 

The  purpose  of  this  thesis  is  to  study  the  thermoacoustic  effect, 
on  a  single  plate  and  on  multiple  plates,  to  better  understand  both  the 
basic  theory  of  thermoacoustic  heat  transport  and  the  reason  for  single 
plate  discrepancies. 

The  results  of  two  experiments  are  reported.  In  the  first  experi¬ 
ment,  changes  in  the  temperature  difference  across  a  TAC  are  studied 
as  a  function  of  the  longitudinal  position  of  the  TAC  in  the  resonant 


8 
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tube,  the  number  of  plates  in  a  stack,  the  spacing  between  plates,  the 
acoustic  pressure  amplitude,  and  the  type  of  gas  in  the  tube.  In  the 
second  experiment,  the  developed  temperature  difference  is  mea¬ 
sured  as  the  TAC  position  is  varied  along  a  diameter  of  the  tube.  The 
results  of  these  experiments  are  compared  to  the  theory  as  developed 
by  Wheatley  and  associates.  Possible  causes  for  the  observed  deviations 
from  theoiy  are  offered. 
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n.  THEORY 


A.  INTRODUCTION 

A  synopsis  of  the  theory  of  thermoacoustic  heat  transport,  based 
upon  John  Wheatley’s  development,  follows.  In  addition  to  the  basic 
theory,  modifications  to  account  for  the  thermal  conductivity  of  the 
gas  will  be  considered.  The  possible  effects  that  non-thermal  pro¬ 
cesses  may  have  on  the  heat  transport  will  also  be  discussed.  The 
reader  is  directed  to  References  1  and  2  for  a  more  detailed  develop¬ 
ment  of  the  basic  theory. 

B.  THERMOACOUSTIC  EFFECT 

1.  Basic  Theory 

In  order  to  understand  the  thermoacoustic  effect,  it  is  neces¬ 
sary  to  describe  what  happens  to  a  gas  parcel  in  an  acoustic  standing 
wave  developed  within  a  resonant  tube.  In  this  study,  we  are  consid¬ 
ering  the  heat  pumping  action  of  such  an  acoustic  wave. 

The  four  stages  of  the  heat  pumping  cycle  are  depicted  in 
Figure  1  [Ref.  2:p.  13].  Initially,  let  us  assume  that  the  plate,  or  TAC, 
has  a  small  temperature  gradient,  VT,  across  it.  Also,  let  us  assume 
that  the  idealized  plate  does  not  conduct  heat  in  the  longitudinal  x- 
direction.  At  step  1,  a  gas  parcel  undergoes  an  adiabatic  compression 
from  pressure  p  to  p+pi.  At  the  same  time,  the  parcel  is  displaced 
toward  the  pressure  antinode  by  a  distance  xj.  This  distance  may  be 
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expressed  as  xi  =pi/pmcw  where  pi  is  the  acoustic  pressure  or  adia¬ 
batic  pressure  change,  pm  is  the  density  of  the  gas,  c  is  the  speed  of 
sound  in  the  gas,  and  co  is  the  angular  frequency.  Also,  at  this  time, 
the  parcel’s  temperature  increases  from  T  to  T+Ti.  Ti  is  the  adiabatic 
temperature  change  and  will  be  discussed  in  more  detail  in  the  fol¬ 
lowing  pages.  At  step  2,  there  is  a  heat  flow,  Q,  in  the  transverse 
direction,  from  the  hotter  gas  volume  to  the  cooler  plate.  This  heat 
flow  also  results  in  a  temperature  reduction  8T  in  the  gas  volume.  On 
the  return  leg,  step  3,  the  gas  parcel  expands  adiabatically  and  moves 
back  a  distance  xi  but  is  now  at  a  lower  temperature  T-5T.  Therefore, 
at  step  4,  there  is  another  heat  exchange,  Q'f  but  this  time  from  the 
hotter  plate  to  the  cooler  gas.  This  heat  exchange  brings  the  gas  back 
to  the  original  temperature  T  while,  at  the  same  time,  it  reduces  the 
temperature  of  this  end  of  the  plate.  This  heat  pumping  cycle  repeats 
with  the  frequency  of  the  acoustic  wave. 

In  describing  the  heat  flow  in  short  stacks  associated  with 
this  process,  Wheatley,  et  al.,  estimate  the  hydrodynamically  trans¬ 
ported  heat  Q"  as 

Q  —  n  8k  u  i  Cp  pm  5T  ( 1 ) 

[Ref  2:p.  14],  where 

n  =  overall  perimeter  of  the  plate  (measured  transverse  to  the 
flow); 

SK  =  thermal  penetration  depth; 

ui  =  acoustic  velocity  magnitude; 

Cp  =  specific  heat  at  constant  pressure  of  the  gas;  and 


5T  =  change  in  temperature  in  the  gas  due  to  the  heat  flow  process. 
From  Figure  1.  step  2,  it  can  be  seen  that 


5T  =  Ti  -  xi  VT  =  Ti  1  - 


T!/Xl 


Ti.  as  earlier  stated,  is  the  adiabatic  temperature  change,  and  may  be 
expressed  by 

Ti  =  ^pi  (3) 

where 

Tm  =  mean  absolute  temperature;  and 
P  =  isobaric  expansion  coefficient. 

For  an  ideal  gas.  the  product  Tmp  equals  unity. 

The  value  for  VT  which  makes  5T  =  0.  in  equation  2,  is 

known  as  the  critical  gradient  so  that 

VTcrtt=^  (4) 

[Ref.  2:p.  14].  This  temperature  gradient  is  the  largest  that  can  be 
acoustically  induced  on  this  plate. 

Using  Wheatley's  definition  T  s  VT/VTCrit  [Ref.  2:p.  14]  and 
rearranging  the  terms  in  equation  (1),  we  can  express  the  hydrody¬ 
namic  heat  flow  as 


9"  ~  n8K  (Tmp)  pi  U]  (1  -D. 
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For  the  assumed  ideal  plate,  we  will  remain  in  the  heat 
pumping  regime  while  the  heat  flow  is  towards  the  nearest  pressure 
antinode.  This  will  occur  as  long  as  VT  <  VTcrlt  or  T  <  1.  If  VT  =  VTcrlt, 
then  there  is  no  heat  flow. 

Now,  let  us  consider  that  the  TAC  is  not  ideal  and,  thus,  does 
conduct  heat.  There  now  exists  a  return  heat  flow  Qr,  due  to  diffusive 
thermal  conduction  through  the  plate  material.  This  heat  flow  may  be 
expressed  as  Qr  =  k  A  VT,  where  k  is  the  thermal  conductivity  in  the 
direction  of  the  plate  length,  A  is  the  area  of  the  plate  transverse  to 
the  flow,  and  VT  is  the  temperature  gradient  across  the  plate.  The 
net  heat  flow  of  this  system,  Qt,  may  be  expressed  as 

Qt  =  Q"-Qr.  (7) 

or 

Qt  =  [n  8k  (Tm3)  piui  (1  -  D]  -  [k  A  VT).  (8) 

Considering  this  heat  flow,  the  heat  pumping  regime  is  main¬ 
tained  as  long  as  Q"  is  directed  toward  the  nearest  pressure  antinode. 
At  Q"  =  Qr.  we  have  reached  steady  state  and  Qt  =  0.  We  have,  then, 
the  maximum  temperature  gradient,  VT,  of  this  system.  For  a  given 
transverse  area  A,  a  non-zero  k  will  limit  the  maximum  value  of  the 
temperature  gradient  to  that  achieved  when  Q"  =  Qr. 

Now,  consider  the  term  n  8<  (Tmp)  as  constant  for  any  posi¬ 
tion  of  the  TAC  in  the  standing  wave.  This  results  in  Q"  being  depen¬ 
dent  on  the  piui  (1  -  T)  product.  [For  small  acoustic  amplitudes.  VT 


is  small,  r  *  0,  and  thus  ( 1  —  T)  — >  1  and  is  inconsequential.  At  large 
amplitudes.  VT  approaches  VTcrit  and  (1  -  r)  ->  0.  This  will  affect  Q". 
These  two  cases  will  be  discussed  in  more  detail  later  in  this  paper.] 
By  studying  equation  8,  we  can  see  that,  for  a  given  k  and  plate  geom¬ 
etry,  the  temperature  gradient  that  exists  across  the  TAC,  in  the  heat 
pumping  case,  is  dependent  on  the  piui  product. 

The  acoustic  amplitude  pi  and  particle  velocity  ui  vary  with 
position  in  the  standing  wave.  They  may  be  expressed  as 

pi  =  Pq  cos  (kx)  (9) 


u  i  =  Uq  sin  (kx). 


(10) 


where  k  is  the  wave  number,  x  is  the  displacement  from  the  closed 
end  of  the  tube,  and  Po  and  Uo  are  the  magnitudes  of  the  pressure  and 
velocity  at  their  respective  antinodes.  Thus  the  product  piui  may  be 


expressed  as 


(Po  U0) 


sin  (2kx) 


The  temperature  gradient,  therefore,  can  be  expected  to  vary 
sinusoidally.  Figure  2,  although  not  to  scale,  shows  the  relationship 
among  pi,  ui,  and  VT  plus  the  variation  of  VTcnt.  As  is  evident,  the 
temperature  gradient  swings  from  positive  to  negative  values  with 
zeroes  at  either  the  pressure  or  velocity  nodes.  This  brings  us  to  the 
point  that  the  heat  flow  is  towards  the  nearest  pressure  antinode. 
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Dependence  of  Heat  Transfer  Parameters 


This  is  expected  because,  as  a  parcel  of  gas  is  moved,  in  the  positive 
x-direction,  toward  a  pressure  antinode  (Figure  1),  the  parcel  is  com- 
pressionally  warmed  and  thus  will  transfer  heat  to  the  plate.  If  the 
parcel  is  moved  toward  a  pressure  node,  it  experiences  cooling  by 
expansion  and  takes  heat  from  the  plate.  Furthermore,  from  equation 
(8)  and  Figure  2,  we  deduce  that  as  VTcrlt  ->  0  the  temperature  gradi¬ 
ent  itself  must  approach  zero. 

In  our  investigation,  the  plate  and  gas  will  initially  be  at  the 
same  temperature.  At  the  onset  of  the  acoustic  pressure  r  =  0  and. 
therefore,  as  seen  by  equation  (5),  the  heat  flow  is  initially  a  maximum. 
This  will  result  in  a  temperature  gradient  across  the  plate  which  will 
increase  toward  VTcrit.  However,  as  discussed  in  equations  (7)  and 
(8),  VT  will  increase  only  until  steady  state  is  achieved,  i.e.,  Qr  =  Q". 

In  Ref.  1  (pp.  155-157),  Wheatley,  et  al.,  discuss  the  afore¬ 
mentioned  heat  flow  in  a  slightly  different  manner.  There,  the  devel¬ 
opment  is  continued  in  order  to  derive  an  expression  with  which  to 
calculate  the  temperature  difference  developed  across  a  TAC.  That 
expression  is: 

'  (.25P§5K(1 +0^1) 

(pmc  (B12/U  (1  ♦  a))  sln  ,2“x/c) 

AT=i  ,i2) 

(4(Kd2/UpmLTmw  (Y  -  1)  (1  -  <J2))  ‘‘  cosl2rox/c)] 

(Ref.  l:p.  157,  Eqn.  17J,  where 

AT  =  temperature  difference  across  the  plate 

Po  =  peak  acoustic  pressure  of  the  driver 


5k  =  thermal  penetration  depth 
a  =  Prandtl  number  of  the  gas 
Y  =  the  ratio  Cp/Cv  of  the  gas 
L  =  length  of  the  plate 
Tm  =  mean  ambient  temperature 
pm  =  mean  density  of  the  gas 
c  =  speed  of  sound  in  the  gas 
to  =  angular  frequency 

x  =  mean  position  of  TAC  center  measured  from  the  closed  end 
of  the  tube 

k  =  thermal  conductivity  of  the  plate 
d2  =  thickness  of  the  plate. 

It  is  this  equation  and  predictions  based  on  it  that  will  be 
examined  in  this  thesis. 

2.  Surrounding  Gas  Considerations 

In  the  preceding  equation,  the  product  xd^  ignores  the 
thermal  return  path  via  the  gas  surrounding  the  TAC.  This  return  path 
is  generally  the  gas  between  two  plates.  Therefore,  the  thickness  of 
the  gas  path  is  equal  to  the  plate  separation.  With  a  stack  of  plates, 
this  path  is  well  defined  and  can  be  expressed  as  either  the  absolute 
distance  between  the  plates  or  as  a  multiple  of  the  thermal  penetra¬ 
tion  depth.  Thus,  we  can  express  the  effective  xd2  for  a  stack  of  plates 


as: 


xd2  =  (Kg  *  separation)  +  (kp  *  d2). 


(13a) 


or 


Kd2  =  (Kg  *  8K  *  N)  +  (kp  *  d2), 


(13b) 
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where  N  is  the  number  of  thermal  penetration  depths  equivalent  to 
the  required  separation  and  Kg  and  tcp  are  the  thermal  conductivities  of 
the  gas  and  plate,  respectively. 

For  a  single  plate,  the  return  path  is  not  as  well  defined.  In 
order  to  estimate  the  effect  of  the  thermal  conduction  through  the 
gas.  we  consider  two  extreme  cases.  In  the  first  case,  the  return  path 
is  assumed  to  be  5  thermal  penetration  depths  above  and  below  the 
plate,  for  a  total  of  10  penetration  depths.  In  the  second  case,  it  is 
assumed  ..nat  the  nearest  boundary  is  approximately  one  tube  radius 
away  and  thus  the  effective  thickness  of  the  return  path  is  one  tube 
diameter.  In  other  words,  the  majority  of  the  space  in  the  tube  pro¬ 
vides  the  return  path.  Using  this  assumption,  the  effective  xd2 
becomes 

xd2  =  (xg  *  Diameter)  +  (kp  *  d2).  (14) 

To  avoid  confusion  between  the  two  methods  of  including  the  gas 
return  path,  equations  (13a)  and  (13b)  will  be  referred  to  as  Theory, 
while  equation  (14)  will  be  referred  to  as  Modified  or  Mod  Theory.  It 
is  predicted  that  Theory  underestimates  the  effect  of  the  gas  thermal 
conduction,  while  Mod  Theory  presents  the  absolute  upper  limit  on 
the  gas  conduction. 

C.  TEMPERATURE  DIFFERENCE  PREDICTION 

Now  consider  a  single  steel  plate  TAC,  of  length  2.26  cm,  width 
2.16  cm,  and  thickness  0.127  mm,  placed  in  an  Argon-filled  tube  at 


294  K.  The  thermal  conductivity  of  steel  is  17  W/m-K  [Ref.  3:p. 
D-187],  The  tube  is  driven  at  705  Hz  at  an  acoustic  pressure  level  of 
150  dB  re  20|iPa,  or  approximately  600  Pa  (rms).  The  gas  return  path 
thickness  is  assumed  to  be  5  thermal  penetration  depths  for  Theory 
and  6  cm  for  Mod  Theory.  Using  these  parameters  and  the  properties 
of  Argon  [Ref.  4;  see  Appendix  11.  the  predicted  temperature  differ¬ 
ence  can  be  calculated  from  equation  (12).  As  can  be  seen  in  Figure  3, 
Theory  predicts  a  higher  temperature  difference  than  Modified 
Theory.  This  difference  is  expected  since  Mod  Theory  assumes  a 
much  larger  volume  of  gas  with  which  to  conduct  heat. 

Consider  now  a  stack  of  plates,  with  plate  separation  of  1  mm, 
placed  in  the  same  resonant  tube  filled  with  Argon  at  294  K.  This 
time,  the  acoustic  pressures  will  be  150,  162,  168,  and  174  dB  re 
20|iPa.  The  predicted  temperature  differences  are  shown  in  Figure  4. 
As  can  be  seen,  the  temperature  difference  developed  across  the  TAC 
increases  with  increasing  acoustic  amplitude.  At  higher  amplitudes, 
there  is  a  tilting  or  skewing  of  the  temperature  curve  toward  the 
pressure  antinode.  In  contrast.  Figure  5  shows  the  effect  of  using 
Helium  in  the  tube.  The  physical  properties  of  this  gas  [Ref.  5:pp  48- 
49]  result  in  a  much  larger  temperature  difference  at  the  same  ampli¬ 
tudes  and  also  lead  to  a  much  more  noticeable  skewing  effect. 

Returning  to  the  discussion  on  heat  flow,  one  can  determine  the 
origin  of  this  skewing.  In  equation  (8)  and  Figure  2,  the  effect  of  the 
product  piui  was  seen.  Near  a  pressure  antinode,  the  pressure  pi  is 
large  and  thus  gas  parcels  undergo  a  large  temperature  change  due  to 


igure  3  -Theory  and  Mod  Theory-Predicted  TAC  Temperature  Difference 


Figure  4  -Steel  Stack  in  Argon;  Spacing=  10  Penetration  Depths 


compression.  At  the  same  time,  the  displacement  due  to  ui  is  small. 
This  combination  of  a  large  temperature  change  over  a  small  dis¬ 
placement  means,  by  definition,  that  there  is  a  large  critical  tempera¬ 
ture  gradient.  This  is  confirmed  by  the  large  value  of  VTcrit  near  the 
pressure  antinodes  in  Figure  2.  This  also  explains  the  large  tempera¬ 
ture  difference  near  the  pressure  antinodes  observed  in  Figures  4  and 
5.  A  large  VTcrit  would  allow  VT  to  increase  via  the  (1  -  T)  term  in 
equation  (8).  and  thus  for  a  given  TAC  length  the  temperature  differ¬ 
ence  would  increase.  Now.  as  the  TAC  is  moved  from  the  pressure 
antinode,  the  value  of  VTcrit  decreases  rapidly.  This  causes  a  rapid 
decrease  in  Q",  because  of  (1  -  T)  ->  0,  and  a  corresponding  rapid 
change  in  the  developed  temperature  gradient  is  predicted.  As  a 
pressure  node  is  approached.  VTcrit  approaches  zero,  but  at  a  slow 
rate.  Therefore,  a  much  more  gradual  change  in  the  temperature  gra¬ 
dient  is  predicted.  This  means  a  corresponding  gradual  change  in  the 
temperature  difference.  Thus,  we  see  that  at  higher  pressure 
amplitudes  we  expect  to  see  the  skewing  evident  in  Figures  4  and  5. 

By  considering  the  separation  of  the  plates  as  a  means  to  define  a 
return  path  through  the  gas.  it  is  evident  that  the  predicted  tempera¬ 
ture  difference  will  vary  with  plate  separation.  Furthermore,  while  it 
is  not  readily  apparent  from  equation  12.  Wheatley’s  group  studied  the 
rate  of  change  in  temperature  on  a  stack  of  plates  as  it  relates  to  plate 
separation  distance  d.  compared  to  the  viscous  penetration  depth 
distance  5y.  Their  conclusion  was  that.  “...  qualitative  changes  from 
the  simple  boundary  layer  theory  may  begin  to  show  up  for  d  <  25v." 


[Ref.  l:p.  155]  We  would,  therefore,  expect  deviation  from  theory  as 
the  separation  approaches  the  viscous  penetration  depth. 

D.  NON-THERMAL  PROCESSES 

In  addition  to  the  aforementioned  thermal  effects,  there  are  other 
hydrodynamic  effects  due  to  the  impinging  acoustic  wave.  While  these 
additional  processes  may  have  secondary  thermal  effects,  they  are 
non-thermal  in  origin  and  therefore  are  identified  as  non-thermal 
processes.  In  a  standing  wave  tube  which  is  driven  at  resonance, 
there  is  a  streaming  motion  at  the  tube  boundaries  which  is  toward 
the  velocity  nodes  as  is  depicted  in  Figure  6  [Ref.  6:p.  348].  As  dis¬ 
cussed  earlier,  the  heat  flow  developed  in  this  same  tube  is  towards 
the  same  velocity  nodes.  If  the  TAC  is  on  the  longitudinal  axis  of  the 
resonant  tube,  the  acoustic  streaming  will  provide  an  increase  in  the 
return  heat  flow.  For  a  single  plate,  this  additional  return  flow  may 
cause  the  temperature  gradient  to  be  less. 

For  a  stack  of  plates,  the  effects  of  streaming  will  depend  on  the 
effect  that  the  additional  plates  have  on  the  streaming  pattern  in  the 
tube.  The  plate  separation  may  also  be  important  to  the  tube  pattern. 
The  one-way  streaming  that  occurs  over  the  individual  plates  may  be 
more  prominent  than  the  tube  effects.  At  this  point,  we  can  only  say 
that  acoustic  streaming  effects  are  difficult  to  predict  but  they  are 
likely  to  occur  at  the  pressure  amplitudes  typically  used  in  the  investi¬ 
gation  of  the  thermoacoustic  effect. 


III.  EXPERIMENTS 


A.  LONGITUDINAL  MAPPING 

In  order  to  study  the  temperature  variation  with  respect  to  posi¬ 
tion  in  a  standing  wave,  a  copper  tube  89.7  cm  in  length  and  6.22  cm 
in  diameter  was  used  as  the  resonant  cavity.  The  tube  was  filled  with 
Argon  and  pressurized  to  just  above  atmospheric  pressure  to  ensure 
that  all  leaks  would  be  outward.  Figure  7  shows  the  equipment 
arrangement  for  this  portion  of  the  experimentation. 

Before  data  acquisition  began,  the  endplate  and  probe  micro¬ 
phones  were  calibrated.  The  microphone  at  the  endplate  (Endevco, 
Model  #8550M1)  was  used  to  accurately  monitor  the  frequency  and 
amplitude  of  sound  in  the  tube.  This  microphone  is  readily  calibrated 
using  a  D.C.  calibration  method  (see  Figure  8).  The  response  of  the 
microphone  was  assumed  to  be  flat  from  D.C.  to  the  experimental  fre¬ 
quencies  because  the  fundamental  microphone  resonance  is  at  50 
KHz.  The  other  microphone  used  was  an  earpiece  from  a  portable 
television  set.  As  this  transducer  was  used  only  to  monitor  the  posi¬ 
tion  of  the  pressure  extremes  in  the  tube  and  not  to  determine  abso¬ 
lute  pressures,  precise  calibration  was  not  required.  Therefore,  a 
rough  comparison  calibration  using  the  Endevco  as  a  standard  was 
completed.  Later,  due  to  poor  performance  at  high  acoustic  pres¬ 
sures.  the  earpiece  was  replaced  by  a  similarly  calibrated,  locally  con¬ 
structed,  piezoceramic  transducer. 
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Figure  B  -  DC  Calibration  of  Endevco  Microphone 


To  determine  a  suitable  driving  frequency,  the  resonances  of  the 
Argon-filled  tube  were  studied  (see  Figures  9  and  10).  An  interesting 
point  of  this  procedure  was  the  discovery  that  the  speed  of  sound  in 
Argon  of  319  m/s  at  0°  C,  obtained  from  a  well-used  source  (Ref.  3:p. 
E-47],  is  in  error.  By  comparing  the  tube  length  to  the  observed  res¬ 
onant  frequencies,  the  speed  of  sound  was  calculated  to  be  near  320 
m/s  at  20°  C.  Further  research  found  three  more-recent  sources  (Ref. 
4;  Ref.  7:p.  3133;  Ref.  8:p.  264]  which  give  the  value  of  319  m/s  at  20° 
C.  This  is  the  value  used  for  the  following  experiments. 

Using  a  single  sheet  glass  TAC,  trial  runs  were  conducted  to  verify 
experimental  procedures.  This  trial  provided  a  method  of  determin¬ 
ing  the  time  required  for  stabilization  of  the  TAC  temperature  after  it 
had  been  moved  (see  Figure  11  for  typical  plot).  This  was  required  for 
the  automation  of  data  acquisition.  Using  the  measured  temperature 
gradient  results  from  this  trial  run  (Figure  12),  a  decision  was  made 
regarding  the  portion  of  the  tube  from  which  useful  data  could  be 
obtained. 

A  TAC,  consisting  of  one  0.0127  cm  thick  stainless  steel  sheet 
laminated  to  one  0.0127  cm  thick  G-10  fiberglass  sheet,  2.26  cm  long 
and  2.16  cm  wide,  was  fitted  with  a  thermopile  consisting  of  four  pairs 
of  constantan-chromel  thermocouples.  To  ensure  correct  voltage 
readings  due  to  the  temperature  difference  developed  across  the  TAC, 
the  leads  from  the  TAC  are  put  in  thermal  contact  with  a  copper  plate 
to  provide  a  stable  temperature  reference  (see  Figure  7).  This  “heat 
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sink"  arrangement  also  limited  random  voltages  caused  by  tempera¬ 
ture  variations  along  the  leads,  but  not  at  the  TAC. 


The  temperature  difference  developed  across  the  TAC,  as  a  func¬ 
tion  of  position  from  the  rigid  end  of  the  tube,  was  measured  for 
acoustic  pressure  levels  of  145,  150,  and  162  dB  re  20  \i  Pa  (rms). 
Next,  stacks  of  3  and  5  plates  were  constructed  by  stacking  one  or  two 
steel  plates  above  and  below  the  measurement  TAC  at  a  separation 
distance  of  approximately  10  thermal  penetration  depths.  Again,  the 
temperature  difference  as  a  function  of  position  was  measured  at 
pressure  levels  150  and  162  dB.  The  measurements  were  then 
repeated  for  a  separation  distance  of  36  thermal  penetration  depths 
and  then  for  a  stack  of  3  plates  separated  by  4  thermal  penetration 
depths. 

To  ensure  that  the  data  were  reproducible,  the  initial  data  mea¬ 
surements  were  taken  twice.  Typical  agreement  in  measured  values  is 
depicted  in  Figure  13.  As  more  data  runs  were  completed,  the  repro¬ 
ducibility  of  the  data  was  verified.  Therefore,  data  from  later  runs 
were  verified  by  retaking  only  random  data  points. 

The  final  measurements  in  this  phase  of  the  study  were  made  in  a 
Helium-filled  tube.  The  new  resonances  of  the  tube  were  found  (see 
Figure  14)  and  the  driving  frequency  2120  Hz  was  chosen  to  maintain 
the  wavelength  to  TAC  length  ratio  used  in  the  Argon-filled  tube.  In 
this  portion  of  the  experiment,  a  stack  of  3  plates  separated  by  6 
thermal  penetration  depths  was  used  for  pressure  levels  of  150  and 
162  dB. 


B.  DIAMETRICAL  MAPPING 


Before  probing  the  tube  diameter,  the  results  of  phase  1  were 
perused  to  determine  a  position  within  the  tube  where  a  reasonable 
temperature  difference  was  established  at  all  pressure  levels  studied. 
This  position  was  41.5  cm  from  the  endplate,  where  2kx  =  11.5.  The 
tube  was  then  cut  at  this  position,  the  apparatus  shown  in  Figure  15 
inserted,  and  the  tube  sections  rejoined.  Again,  the  tube  was  filled 
with  Argon  and  the  resonances  were  found  (see  Figure  16).  The  reso¬ 
nance  frequency  closest  to  the  driving  frequency  used  in  the  first 
experiment  was  chosen  as  the  driving  frequency  for  this  experiment. 
For  acoustic  pressure  levels  of  150  and  162  dB,  the  temperature  dif¬ 
ference  across  the  single  TAC  was  measured  as  a  function  of  position 
along  the  diameter  of  the  tube.  Then  the  same  measurement  was 
made  for  the  center  plate  of  a  stack  of  3  plates  separated  by  10  ther¬ 
mal  penetration  depths. 
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-Insert  to  Measure  Diametric  Temperature  Gradient  Variations 


igure  16-Natural  Resonances  of  Argon  Filled  Tube  with  Insert 


IV.  RESULTS 


A.  LONGITUDINAL  MAPPING 

As  previously  described,  a  TAC  constructed  of  a  single  glass  plate 
was  used  in  a  trial  measurement.  As  can  be  seen  in  Figure  17,  at  an 
acoustic  amplitude  near  150  dB  re  20|iPa  (rms),  the  measured  tem¬ 
perature  difference  and  acoustic  amplitude  have  the  predicted  rela¬ 
tionship  (i.e.,  the  temperature  difference  varies  approximately 
sinusoidally  with  zero  crossings  corresponding  to  acoustic  pressure 
maxima  and  minima). 

The  temperature  difference  across  the  single  stainless  steel  TAC 
was  measured  as  a  function  of  tube  position  and  pressure  amplitudes. 
The  results  are  shown  in  Figures  18  through  20.  In  these  figures,  the 
solid  line  represents  the  experimental  values;  the  dashed  and  dotted 
lines  represent  Theory  and  Mod  Theory,  respectively.  At  acoustic 
amplitudes  of  145  and  150  dB  re  20pPa  (rms).  the  measured  temper¬ 
ature  difference  falls  between  the  values  predicted  by  both  Theory  and 
Mod  Theory  calculation  methods  (Figure  18  and  19).  At  the  higher 
level  of  162  dB.  the  peak  temperature  difference  is  at  least  2°  C  (or 
approximately  50  percent)  lower  than  that  predicted  by  either 
theoretical  method  (Figure  20).  However,  there  is  still  good 
agreement  near  the  pressure  antinode  (2kx  *  12.7).  Additionally, 
skewing  of  the  data  toward  the  pressure  antinodes  is  evident,  yet  not 
predicted. 


When  stacks  of  3  and  5  plates,  with  a  plate  separation  of  36  ther¬ 
mal  penetration  depths,  are  used  at  an  acoustic  pressure  of  150  dB. 
there  again  is  close  agreement  between  measured  and  predicted  data 
(Figure  21).  The  agreement  between  the  observed  temperature  dif¬ 
ferences  for  a  3  or  5  plate  stack  is  within  experimental  fluctuation 
observed  on  repeat  data  runs,  as  discussed  earlier  (Figure  13).  At  the 
higher  sound  level  of  162  dB,  a  large  discrepancy,  2°  C  (about  50  per¬ 
cent),  between  measured  and  predicted  peak  temperature  differences 
is  observed  (Figure  22).  although  the  agreement  near  the  pressure 
antinode  is  still  good.  Again,  the  3  and  5  plate  stacks  produce  compa¬ 
rable  results. 

When  measurements  are  made  using  stacks  with  plate  separations 
of  10  thermal  penetration  depths,  there  is  again  agreement  between 
measured  and  predicted  data  at  the  low  acoustic  amplitude  (Figure  23) 
and  disagreement  at  the  high  amplitude  (Figure  24).  In  this  case,  the 
magnitude  of  the  disagreement  is  smaller  (only  1°  C  or  about  30  per¬ 
cent)  than  the  predicted  peak.  This  disagreement  is  not  as  much  as 
that  for  a  single  plate  or  for  a  stack  with  36  penetration  depth  spacing. 
Additionally,  there  again  appears  to  be  a  strong  similarity  between  the 
results  for  the  3  plate  and  5  plate  stacks. 

With  a  stack  consisting  of  3  plates,  separated  by  4  thermal  pene¬ 
tration  depths,  again  there  is  good  agreement  between  observed  and 
theoretically  calculated  temperature  differences  at  a  low  level  of 
acoustic  pressure  (Figure  25).  Figure  26  shows  the  now-familiar  dis¬ 
crepancy  at  the  higher  driving  amplitude,  but  it  has  increased  back  to 


Figure  23  -Steel  Stack;  Spacing=iO  Penetration  Depths;  SL=150dB 
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Stack;  Spacing*5  4  Penetration  Depths;  SL=162dB 


approximately  2°  C  (about  50  percent).  The  noticeable  skewing  is 
again  present,  with  good  agreement  near  the  pressure  antinode. 

Finally,  this  same  TAC  was  used  in  a  Helium-filled  tube  with  a 
plate  separation  of  6  thermal  penetration  depths.  While  the  measured 
data  are  approximately  1°  C  (50  percent)  lower  than  predicted  at  the 
150  dB  level,  the  shape  of  the  resulting  curve  is  as  predicted  (Figure 
27).  However,  at  the  162  dB  level,  the  discrepancy  increases  to 
approximately  16°  C  (about  70  percent)  lower  than  the  maximum  pre¬ 
dicted.  Also,  as  with  the  Argon  data,  the  agreement  is  best  near  the 
pressure  antinode.  Furthermore,  the  observed  skewing  is  much  more 
pronounced  than  predicted  (Figure  28).  In  this  figure,  the  encircled 
jump  in  the  data  is  the  result  of  a  blown  fuse  in  the  driver  circuit 
which  had  to  be  replaced  before  data  acquisition  could  continue. 

B.  DIAMETRICAL  MAPPING 

In  this  experiment,  the  diameter  of  the  tube  was  probed  with  two 
TAC  structures.  First,  using  a  single  steel  TAC.  the  temperature  dif¬ 
ference  variations  across  the  diameter  were  measured  for  acoustic 
amplitudes  of  150  and  162  dB  (Figures  29  and  30).  At  the  lower  level, 
the  temperature  difference  magnitude  increases,  albeit  less  than 
0.01°  C  (3  percent),  as  the  TAC  was  moved  across  the  diameter.  At  the 
162  dB  level,  the  difference  magnitude  increases  about  0.15°  C  (10 
percent)  across  the  diameter.  The  probe  construction  prohibited 
measurement  within  0.9  cm  of  the  tube  wall.  The  3.11  cm  distance 
corresponds  to  the  tube  center. 


Figure  28  -Steel  Stack  in  Helium;  Spacingc6  Penetration  Depths;  SL-162dB 
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Figure  29-  Diametric  Temperature  Variation-Single  TAC-150dB 


Figure  30  -  Diametric  Temperature  Variation-Single  TAC-162dB 


Next,  the  measurements  were  repeated  using  a  3  plate  stack  with 
a  plate  separation  of  10  thermal  penetration  depths  (Figures  31  and 
32).  At  150  dB,  there  is  a  slight  decrease  in  magnitude,  0.05°  C  (about 
10  percent),  across  the  diameter.  There  was  an  interesting  result 
observed  at  162  dB:  The  temperature  difference  went  through  a 
maximum  near  the  tube  center,  with  a  temperature  swing  of  0.2°  C 
(5  percent). 


Figure  32  -  Diametric  Temperature  Variation-Steel  Stack-162dB 


V. 


The  modifications  made  to  Wheatley’s  theory  allow  for  considera¬ 
tion  of  the  return  path  through  the  gas,  which  Wheatley  considered 
insignificant  for  a  steel  plate  TAC.  For  a  single  plate  TAC  at  low  ampli¬ 
tude.  Theory  and  Mod  Theory  are  both  in  good  agreement  with  the 
measured  data  (Figures  18  and  19).  It  was  stated  in  Chapter  II  that 
the  Theory  and  Mod  Theory  represented  lower  and  upper  limits  to 
the  temperature  difference  across  the  TAC.  This  prediction  was  con¬ 
firmed  at  low  pressure  amplitudes  where  the  experimental  data  falls 
between  the  two  predictions  (Figures  18  and  19).  At  the  higher 
acoustic  amplitude,  neither  theory  can  account  for  the  lower  tempera¬ 
ture  gradient  observed  across  the  TAC  (Figure  20). 

It  was  believed  that  by  studying  stacks  of  plates  the  cause  of  this 
discrepancy  could  be  found.  For  stacks  of  3  and  5  plates  at  various 
separations,  the  measured  data  at  low  acoustic  amplitudes  compares 
favorably  with  theory  (Figures  21.  23,  and  25).  This  agreement  is 
observ  ed  no  matter  the  spacing  or  the  number  of  plates.  Even  with  a 
different  gas.  there  is  reasonable  agreement  (Figure  27).  However,  at 
the  higher  acoustic  amplitude,  the  magnitude  of  the  temperature  dif¬ 
ference  is  noticeably  smaller  than  predicted,  regardless  of  the  separa¬ 
tion.  the  number  of  plates,  or  the  gas  in  the  tube  (Figures  22,  24,  26, 
and  28).  Additionally,  the  skewing  phenomenon  is  observed.  It  is. 
therefore,  evident  that  the  inclusion  of  return  path  through  the  gas  in 
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the  theory  is  not  sufficient  to  account  for  the  discrepancies  observed. 
Evidently,  at  higher  amplitudes  some  other  processes  are  causing  the 
discrepancies. 

It  is  interesting  to  note  the  effect  that  variations  in  plate  separa¬ 
tion  have  on  the  experimental  results.  As  the  plate  separations  are 
reduced  from  3  cm,  i.e.,  a  single  plate,  to  a  separation  of  10  thermal 
penetration  depths,  the  agreement  with  theory  improves.  But  at  a 
smaller  separation  distance  (i.e.,  4  penetration  depths),  the  agreement 
begins  to  diverge  again  (Figure  33).  As  the  separation  distance 
decreases  from  3  cm  to  4  thermal  penetration  depths,  the  observed 
temperature  difference  amplitude  peak  first  increases  then  decreases. 
Concurrently,  the  slope  of  the  temperature  difference  curve  first 
steepens,  then  becomes  more  shallow.  This  indicates  that  there 
exists  an  optimum  plate  spacing  distance. 

The  observed  skewing  is  predicted  by  the  theory  (Figures  4  and 
5)  but  at  much  higher  amplitudes  than  those  used  in  these  experi¬ 
ments.  In  addition,  theory  predicts  that  the  skewing  is  more  pro¬ 
nounced  in  Helium  than  in  Argon.  Comparison  of  the  results  for  these 
two  gases  at  162  dB  pressure  amplitude  reveals  that  the  skewing  is 
equally  severe.  This  fact,  coupled  with  the  fact  that  the  onset  of 
skewing  occurs  at  much  lower  acoustic  amplitudes  than  predicted, 
indicates  that  the  reasons  for  the  observed  skewing  are  different  from 
those  presented  earlier  in  this  thesis. 

It  should  be  pointed  out  that  the  agreement  between  theory  and 
experiment  is  quite  good  at  all  pressure  amplitudes  in  regions  near 
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the  pressure  antinode.  In  this  region,  the  velocity  amplitude  is  small 
and  r  is  small.  The  discrepancy  is  worst  in  regions  of  large  velocity 
amplitude  or  large  r  (=  1).  This  observation  points  to  possible  regions 
of  further  investigation  in  order  to  explain  the  discrepancy. 

Theory  assumes  that  the  TAC  is  short  compared  to  the  wave¬ 
length,  so  that  the  temperature  gradient  is  constant  along  the  plate 
and  that  acoustic  dissipation  due  to  the  plate  is  small.  In  our  mea¬ 
surements,  the  TAC  was  approximately  5  percent  of  the  wavelength. 
Judging  by  the  close  agreement  between  the  results  and  theory,  the 
length  of  the  TAC  is  sufficiently  short  for  low  amplitudes.  At  higher 
amplitudes,  the  short  stack  assumptions  may  be  invalid,  resulting  in 
the  observed  discrepancies. 

Still  another  possible  reason  for  the  discrepancies  could  be  that 
the  harmonics  of  the  driving  frequency  are  influencing  the  predicted 
heat  flow.  Using  equation  (12)  and  the  acoustic  amplitude  of  the 
second  harmonic,  measured  to  be  some  20  dB  lower  than  that  of  the 
resonance  frequency,  the  temperature  difference  contribution  was 
calculated.  As  can  be  seen  in  Figure  34,  the  maximum  temperature 
difference  due  to  the  harmonic  is  only  0.04°  C.  Therefore,  the  second 
harmonic  is  not  the  cause  of  the  discrepancies. 

If  we  consider  the  theory  as  correct,  and  if  it  includes  all  the 
major  thermoacoustic  considerations,  then  some  processes  of  non- 
thermal  origin  are  responsible  for  the  observed  discrepancies  from 
theory.  One  such  non-thermal  process  is  acoustic  streaming. 


HARMONIC  FREQ  -  1410  Hz 


d  Harmonic  Contribution 


Because  acoustic  streaming  occurs  over  any  surface  affected  by  an 
acoustic  wave,  we  must  consider  two  sources  which  could  influence 
the  heat  transport.  The  first  source  is  the  streaming  associated  with 
the  tube  walls.  This  flow  pattern  is  as  discussed  earlier  (Figure  6). 
The  second  source  is  the  streaming  generated  by  the  TAC  plates 
themselves.  It  is  assumed  that  the  flow  generated  by  the  streaming 
would  transport  heat  convectively,  similar  to  what  would  occur  if  a 
stream  of  gas  were  blown  across  a  plate  having  a  temperature  gradient 
across  it.  If  the  streaming  effects  of  the  tube  are  dominant,  then  the 
TAC  temperature  difference  should  vary  along  a  diameter  of  the  tube. 
This  variation  should  be  manifested  as  a  difference  minimum  at  the 
tube  center  and  at  tube  walls,  with  a  difference  maximum  midway 
between  them.  The  measured  diameter  variations  in  the  temperature 
difference  do  not  demonstrate  such  behavior. 


The  measured  data  for  both  a  single  plate  TAC  and  a  stack,  at  both 
150  and  162  dB  acoustic  pressure  levels,  is  essentially  constant, 
within  experimental  error,  across  the  diameter  (Figures  29-32).  The 
magnitude  of  experimental  error  was  determined  by  comparing  the 
magnitude  of  the  variations  to  those  in  successive  data  runs  (Figure 
13).  and  to  the  stability  of  the  measured  temperature  differences 
(Figure  11).  It  is.  therefore,  concluded  that  acoustic  streaming  gener¬ 
ated  by  the  tube  is  not  the  cause  of  the  discrepancies.  Streaming 
effects  on  the  plates  occur  within  a  few  viscous  penetration  depths. 
The  viscous  penetration  depth  in  these  measurements  is  on  the  order 
of  70  jim.  The  thermocouple  wires  are  3  mils  or  76.2  nm  in  diameter, 
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so  the  thermocouple  joints  are  larger  than  the  viscous  penetration 
depth.  Therefore,  the  effects  of  plate  streaming  could  not  be 
determined. 

The  changes  in  the  observed  data,  when  compared  to  the  previ¬ 
ously  measured  results,  indicate  an  apparent  motion  of  the  TAC  toward 
the  driver.  Even  when  the  TAC  is  purposely  positioned  1  cm  closer  to 
the  driver,  there  still  exists  an  apparent  motion  (Figures  35  and  36). 
The  observed  difference  in  the  slopes  in  Figure  36  is  due  to  the  TAC 
starting  at  positions  on  opposite  sides  of  the  peak  of  the  temperature 
gradient,  which  exists  at  the  probe  position.  To  determine  the  cause 
of  this  apparent  motion,  the  precision  of  the  probe  and  tube  construc¬ 
tion  were  studied,  as  were  the  variations  in  driving  frequency,  driving 
amplitude,  and  ambient  temperature.  These  changes  are  considered 
to  be  of  little  consequence  and  the  exact  reasons  for  the  apparent 
motion  are  undetermined. 
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Diametric  Temperature  Variation-Single  TAC-150dB 


Figure  36  -  Diametric  Temperature  Variation-Single  TAC-162dB 


It  is  concluded  that  the  theory,  for  both  a  single  plate  and  stack 
configurations,  is  valid  at  low  acoustic  amplitudes  but  not  at  higher 
amplitudes.  Furthermore,  assuming  that  all  thermal  processes  have 
been  included  in  the  present  theory,  it  is  concluded  that  non-thermal 
processes,  other  than  acoustic  streaming  in  the  resonant  tube,  are 
responsible  for  deviations  from  the  theory.  Finally,  because  of  the  dis¬ 
crepancies  observed,  especially  at  high  acoustic  amplitudes,  it  is  con¬ 
cluded  that  a  TAC  would  be  a  poor,  or  at  least  a  very  limited,  acoustic 
power  probe,  as  suggested  by  Wheatley  [Ref.  2:p.  121,  until  a  better 
understanding  of  high  amplitude  effects  is  achieved. 

To  better  understand  the  processes  involved  and  to  determine 
the  cause  of  the  discrepancies,  it  is  recommended  that  the  experi¬ 
ments  be  repeated  using: 

a  Helium  to  maximize  the  thermoacoustic  effects: 

b.  a  shorter  TAC  to  give  better  agreement  with  the  short  plate 
theory; 

c.  a  phase-lock  loop  to  stabilize  the  driving  frequency:  and 

d.  various  plate  separations  in  a  stack  to  determine  the  optimum 
plate  separation  as  a  function  of  thermal  penetration  depth. 


APPENDIX 


PARTIAL  LISTING  OF  THE  PHYSICAL  PROPERTIES 
OF  ARGON  AND  HELIUM 
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